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Abstract 

This  report  details  the  release  of  version  2.3  data  for  CHAMP  and  GRACE  accelerometer-derived 
thermospheric  total  mass  densities.  The  relevant  references  and  updated  procedures  with  respect  to 
version  2.2  data  are  cited. 

Description  of  Data 

The  CHAMP  [Reigber  et  al.,  2002]  and  GRACE  [Tapley  et  al.,  2004]  satellites  have  provided  useful  data  to 
the  upper  atmospheric  and  ionospheric  scientific  community  since  2000  and  2002,  respectively.  The 
accelerometers  onboard  these  satellites  measure  the  vector  quantity  of  acceleration  caused  by  non- 
gravitational  forces.  After  modeling  and  removal  of  the  acceleration  signals  caused  by  solar  radiation, 
Earth's  albedo  and  infrared  radiation,  the  drag  acceleration,  aD,  can  be  expressed  in  the  direction  of  the 
satellite  velocity  with  respect  to  the  atmosphere,  v: 

aD  =  ~^pCDAref\v\2  (1) 

where  p  is  the  total  atmospheric  density,  CD  is  the  coefficient  of  drag  and  Aref  is  the  reference  area  of  the 
satellite.  Similarly,  the  lift  acceleration,  aL,  can  be  expressed  in  a  direction  perpendicular  to  v  that 
depends  on  both  the  satellite  geometry  and  the  assumed  drag/lift  force  model: 

aL  = -\pCLAref\v\2  (2) 

where  CL  is  the  coefficient  of  lift.  These  equations  can  be  directly  solved  for  atmospheric  density, 
employing  certain  assumptions: 

1.  The  measured  acceleration  contains  only  the  signals  related  to  atmospheric  drag. 

2.  The  coefficients  of  drag  and  lift  can  be  modeled.  The  current  version  adheres  to  the  method  of 
Sutton  [2009], 

3.  The  reference  area  is  known  from  knowledge  of  the  satellite  geometry  and  attitude. 

4.  Atmospheric  winds  are  known  or  can  be  neglected. 

In  addition,  an  estimate  of  the  error  in  density  caused  by  our  limited  knowledge  of  the  parameters  in 
equations  1-2,  as  well  as  in  the  modeled  accelerations  caused  by  solar  radiation  pressure,  Earth's  albedo 
and  infrared  radiation,  can  readily  be  calculated.  The  uncertainty  in  density,  up,  can  be  expressed  using 
the  following  formula: 
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where  the  x/s  are  the  n  parameters  introducing  error  into  the  measurement  and  the  ux.'s  are  the 
estimated  uncertainties  of  each  parameter.  In  our  analysis,  the  ux.'s  are  assumed  to  be  uncorrelated. 
For  more  details  on  the  data  reduction  and  error  analysis  process  see  Sutton  et  al.,  2007. 

In  an  effort  to  reduce  the  size  of  the  density  data  set  without  degrading  the  quality,  the  data  has  been 
binned  and  averaged  along  the  satellite's  orbit  in  3-degree  increments.  As  the  sample  size  within  a  bin, 
n,  increases,  the  constant  systematic  errors  tend  to  average,  while  random  errors  tend  to  decrease  in 
proportion  to  1/Vri.  This  implies  that  it  is  possible  to  minimize  the  error  of  the  averaged  parameter  by 
discriminating  which  data  to  include  in  the  average.  The  algorithm  responsible  for  this  step 
automatically  finds  the  optimal  combination  of  data  that  minimizes  the  combined  error  of  the  averaging 
bin,  and  discards  the  rest.  In  practice,  however,  only  a  small  percentage  of  density  data  was  actually 
discarded. 

In  the  event  that  binned  and  averaged  data  does  not  meet  the  requirements  of  your  work,  please 
contact  Dr.  Eric  Sutton  at  "AFRL/RVBXI  Org  Mailbox  <afrl.rvbxiorgmailbox(5)kirtland.af.mil>"  to  request 
higher  resolution  data  (0.1  Hz  for  CHAMP  and  0.2-1  Hz  for  GRACE). 

Changes  for  Version  2.3 

The  major  changes  for  this  version  are  summarized  here: 

1.  Extension  of  the  data  sets  for  CHAMP  and  GRACE. 

2.  The  NRLMSIS  model  [Picone  et  al.,  2002]  is  now  run  using  FlO.7  and  Fio.7  values  as  measured  at 
the  surface  of  the  Earth.  Data  prior  to  version  2.3  used  Fw.  7  and  Fw.  7  values  normalized  to  1  AU 
which  further  complicated  any  data-model  comparisons. 

3.  Computation  of  the  quasi-dipole  latitude,  magnetic  longitude,  and  magnetic  local  time  was 
updated  to  the  more  compact  and  efficient  routines  given  by  Emmert  et  al.  [2010], 

Terms  of  Use 

These  data  are  processed  and  made  publicly  accessible  as  a  service  to  the  scientific  community. 
Individuals  planning  scientific  investigations  with  these  data  should  contact  Dr.  Eric  Sutton  at 
"AFRL/RVBXI  Org  Mailbox  <afrl.rvbxiorgmailbox(5)kirtland.af.mil>"  to  verify  the  status  of  possible  new 
versions  of  the  data,  and  to  arrange  coordination  and/or  collaboration  with  our  ongoing  scientific 
activities,  as  appropriate.  Please  reference  the  current  document  in  all  publications  and  presentations: 

Sutton,  E.  K.  (2011),  Accelerometer-Derived  Atmospheric  Densities  from  the  CHAMP  and  GRACE 
Accelerometers:  Version  2.3,  AFRL  Technical  Memo.  DTIC#  ADxxxxxx. 
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Future  Plans 

In  future  data  versions,  we  plan  to  implement  an  updated  method  for  modeling  the  drag  and  lift 

coefficients  to  account  for  shadowing  of  the  neutral  atmosphere  by  the  satellite  geometry  and  multiple 

reflections  of  atmospheric  particles. 
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Appendix  A:  Sample  Plots 

The  following  are  sample  time  series  plots  of  GRACE-B  density  and  the  corresponding  NRLMSIS  model 
density  for  2009  and  the  first  half  of  2010: 
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Figure  1:  GRACE-B  density  (blue)  and  NRLMSIS  (green)  sampled  on  the  satellite's  orbit  for  2009. 
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Figure  2:  Same  as  figure  1,  for  the  first  half  of  2010. 
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Appendix  B:  Example  Ascii  File  Structure 

Version  2.3;  created  by  suttonek  on  30-Nov-2010  21:14:56 

Two-digit  Year  (years) ;Day  of  the  Year  (days ); Second  of  the  Day  (GPS  time, sec) ; Center  Latitude  of  3-degree  Bin 
(deg) ; Satellite  Geodetic  Latitude  (deg) ; Satellite  Longitude  (deg) ; Satellite  Height  (km) ; Satellite  Local  Time 
(hours ); Satellite  Quasi-Dipole  Latitude  (deg) ; Satellite  Magnetic  Longitude  (deg) ; Satellite  Magnetic  Local  Time 
(hours ); Neutral  Density  (kg/mA3) /Neutral  Density  Normalized  to  400km  using  NRLMSISeOO ;  Neutral  Density  Normalized 
to  410km  using  NRLMSISeOO ; NRLMSISeOO  Neutral  Density  at  Satellite  Height; Uncertainty  in  Neutral  Density 
(kg/mA3) /Number  of  Data  Points  in  Current  Averaging  Bin;Number  of  Points  in  Current  Averaging  Bin  that  were 
affected  by  Thrusters; Average  Coefficient  of  Drag  Used  in  Current  Averaging  Bin 


10  177 
10  177 
10  177 
10  177 
10  177 
10  177 
10  177 
10  177 
10  177 
10  177 
10  177 
10  177 
10  177 


10.00 

23.33 

47.50 

85.00 

162.50 
214.00 
250.00 
300.00 

347.50 

392.50 
440.00 

497.50 
535.00 


-90 

-90 

-87 

-84 

-81 

-78 

-75 

-72 

-69 

-66 

-63 

-60 

-57 


-88.99419 

-88.84385 

-87.64724 

-85.38886 

-80.55098 

-77.31195 

-75.04370 

-71.88975 

-68.89044 

-66.04656 

-63.04223 

-59.40206 

-57.02601 


-63.98143  503.413  19.7332  -73.71119 
-18.88273  503.416  22.7435  -73.77020 
14.34620  503.392  0.9654  -73.77390 


26.54422  503.274 
32.56386  502.721 
33.89011  502.131 
34.42544  501.615 
34.88957  500.764 
35.14952  499.815 
35.29419  498.798 
35.37486  497.607 
35.40426  496.019 
35.39422  494.905 


1.7891  -73.52399 
2.2119  -72.14352 
2.3146  -70.72843 
2.3603  -69.57310 
2.4051  -67.80791 
2.4357  -66.01116 
2.4578  -64.23829 
2.4764  -62.31757 
2.4943  -59.94606 
2.5041  -58.37744 


16.49485  20.5873  7.108051^14  5.098546^13  4.105823^13  9.019177^14  3.860230^15 

19.37784  20.7837  7.097737^14  5.097750(^13  4.111139^13  9.498847^14  3.858130^15 

24.59991  21.1394  7.089113^14  5.090842^13  4.104115^13  9.416610^14  3.862704^15 

32.52198  21.6793  7.166710^14  5.142491^13  4.143259^13  9.302305^14  3.937020^15 

47.24453  22.6851  7.210448^14  5.153559^13  4.146200^13  9.102788^14  3.934854^15 

55.38685  23.2441  6.515830^14  4.636213^13  3.725958^13  8.989336^14  3.625738^15 

60.29750  23.5828  6.212690^14  4.402328e-13  3.535183(^13  8.916081^14  3.469878^15 

66.17679  23.9904  5.887052^14  4.141842(^13  3.322204^13  8.820172^14  3.327229^15 

70.91245  0.3211  5.725469(^14  3.994563^13  3.200514(^13  8.733658^14  3.266365^15  8 

74.78696  0.5935  5.472762^14  3.782755(^13  3.027602^13  8.655030^14  3.141633^15  10 

78.36586  0.8471  5.472776(^14  3.740401^13  2.990373^13  8.575618e-14  3.149964^15  9 

82.15424  1.1177  5.220283^14  3.513384^13  2.805144(^13  8.486031^14  3.053279^15  6 

84.37216  1.2774  5.160748(^14  3.435582^13  2.740690e-13  8.432946=5-14  2.998837^15  9 


0  3.567 

1  3.556 

2  3.561 
1  3.569 
0  3.584 
0  3.554 
0  3.528 
0  3.495 
0  3.459 
0  3.425 
0  3.391 
0  3.351 
0  3.331 


-  1843  Lines  Not  Shown  - 


10  177 

85897.50 

-3 

-2.03912 

40.08625 

465.171 

2.5287 

-10.86710 

113.14250 

2.8705 

5.774891(^14 

2.616695^13 

2.059126^13 

1.099009^13 

3.289992^15 

4 

0 

2.865 

10  177 

85934.17 

0 

0.31897 

39.97335 

464.313 

2.5313 

-8.33138 

113.15544 

2.8828 

5.492340e-14 

2.456696^13 

1.931637(^13 

1.099238&-13 

3.146030^15 

6 

0 

2.857 

10  177 

85985.00 

3 

3.58904 

39.81690 

463.224 

2.5350 

-4.79005 

113.10157 

2.8951 

5.272951(^14 

2.321755e-13 

1.823306^13 

1.092631e-13 

3.049638^15 

4 

0 

2.846 

10  177 

86037.50 

6 

6.96718 

39.65576 

462.228 

2.5389 

-1.10935 

112.97285 

2.9029 

5.043097^14 

2.190352(^13 

1.717852(^13 

1.078722e-13 

2.941714^15 

4 

0 

2.840 

10  177 

86075.00 

9 

9.38049 

39.54115 

461.597 

2.5416 

1.52951 

112.84460 

2.9061 

4.742164e-14 

2.042500^13 

1.600410^13 

1.065679=5-13 

2.805641^15 

4 

0 

2.837 

10  177 

86097.50 

12 

10.82857 

39.47267 

461.250 

2.5433 

3.11540 

112.75604 

2.9072 

4.851799e-14 

2.080318^13 

1.629176e-13 

1.057134^13 

2.874003^15 

2 

0 

2.835 

10  177 

86162.50 

15 

15.01202 

39.27634 

460.383 

2.5483 

7.70243 

112.46505 

2.9082 

5.313021(^14 

2.252150e-13 

1.761382e-13 

1.032358&-13 

3.086137^15 

2 

2 

2.835 

10  177 

86225.00 

18 

19.03435 

39.09039 

459.735 

2.5532 

12.11338 

112.15878 

2.9073 

4.966526e-14 

2.084956^13 

1.629274e-13 

1.013548^13 

2.920167^15 

3 

0 

2.839 

10  177 

86242.50 

21 

20.16050 

39.03895 

459.584 

2.5547 

13.34746 

112.07222 

2.9071 

4.928709e-14 

2.063579^13 

1.612379^13 

1.010017e-13 

2.896263^15 

4 

0 

2.841 

10  177 

86309.17 

24 

24.44990 

38.84606 

459.134 

2.5603 

18.04157 

111.75449 

2.9068 

4.281129e-14 

1.773611(^13 

1.385968e-13 

1.006618^13 

2.608541^15 

6 

0 

2.844 

10  177 

86335.00 

27 

26.11168 

38.77285 

459.009 

2.5626 

19.85662 

111.64132 

2.9073 

4.261275e-14 

1.757480e-13 

1.373793^13 

1.010458&-13 

2.633609^15 

3 

0 

2.843 

10  177 

86385.00 

30 

29.32731 

38.63407 

458.840 

2.5673 

23.36192 

111.44787 

2.9101 

5.680612e-14 

2.320233^13 

1.815635^13 

1.027581e-13 

3.357270^15 

2 

2 

2.848 
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Appendix  C:  Support  Files 

In  addition  to  this  document,  two  support  text  files  accompany  the  data.  The  'HISTORY. density.v2. 3.txt' 
file  outlines  the  changes  that  have  been  made  to  the  processing  routines  from  version  to  version.  The 
'README. density. v2. 3.txt'  gives  instructions  for  acquiring  this  data  and  interpreting  the  data  format. 

History  File: 

Changes  for  Version  2.0  CHAMP  Density: 


1.  Interpolation  of  satellite  attitude  data  improved 

a.  methods  updated  to  handle  missing  data 

b.  index  (Numlnterp  variable)  added  to  data  set  indicating  any 
interpolation  of  attitude  data 

2.  Accelerometer  instrument  biases  reprocessed  and  smoothed 

3.  Coefficient  of  drag  updated  from  method  of  Cook  [1965]  to  Sentman  [1961] 

a.  includes  more  realistic  formulation  of  the  coefficient  of  drag 

b.  includes  the  effect  of  thermal  drag,  which  increases  the  coefficient 
of  drag  for  long  satellites 

See:  Sutton,  E.  K.  (2009),  Normalized  Force  Coefficients  for 

Satellites  with  Elongated  Shapes,  J.  Spacecraft  and  Rockets, 
46(1),  112-116,  doi:10. 2514/1. 40940. 

c.  Variable  'Cd'  added 

4.  Variable  'LatBin'  added  and  ' Lat '  modified 

a.  LatBin  gives  the  center  of  the  current  averaging  bin 

b.  Lat  gives  the  mean  position  of  the  satellite  at  data  points  within 
the  current  averaging  bin 

c.  This  ensures  that  the  satellite  position  (defined  by  variables  Lat, 
Lon,  and  Height)  is  consistent 

5.  NRL-MSIS  model  density  is  given  at  each  time  and  satellite  position 


Changes  for  Version  2.1  CHAMP  Density: 


1.  Addition  of  Quasi-Dipole  Latitude  (Mlat) ,  Magnetic  Longitude  (Mlon) ,  and 
Magnetic  Local  Time  (Mlt) 

See:  Richmond,  A.  D.  (1995),  Ionospheric  Electrodynamics  Using 

Magnetic  Apex  Coordinates,  J.  Geomag.  Geoelectr.,  47(2), 
191-212 . 


Changes  for  Version  2.2  CHAMP  and  GRACE  Density: 


1.  First  Version  for  GRACE  Density,  everything  conforms  to  version  2.2  of 
CHAMP  density 


Changes  for  Version  2.3  CHAMP  and  GRACE  Density: 


1.  NRL-MSIS  density  values  updated 

a.  now  using  observed  F10.7  values  in  place  of  1  AU  normalized  values 

b.  corrected  Ap/ap/F10.7  indices  used  during  the  first  couple  data 
points  of  the  day 

2.  Software  update  for  computation  of  Quasi-Dipole  Latitude,  Magnetic 
Longitude,  and  Magnetic  Local  Time 

See:  Emmert,  J.  T.,  A.  D.  Richmond,  and  D.  P.  Drob  (2010),  A 
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computationally  compact  representation  of  Magnetic  Apex  and 
Quasi  Dipole  coordinates  with  smoothe  base  vectors,  J. 
Geophys .  Res.,  115,  A08322,  doi : 10 . 1029/2010JA15326 . 


Readme  File: 

%%%%  Format  Description: 

All  Density  data  has  been  averaged  into  3-degree  latitudinal  bins  to 
conserve  space  and  reduce  any  random  errors.  The  bins  are  centered  around  lat 
=  [-90:3:90]  deg,  omitting  any  values  for  which  there  are  0  data  points  in 
the  averaging  bin.  Neutral  density  has  also  been  normalized  to  2  heights 
(400km  and  410km)  using  the  NRL-MSISe-00  Empirical  Density  Model. 


%%%%  Ascii  Format: 

Ascii  files  are  arranged  in  columnn  format  with  a  2-line  header  of 
Version  Information  and  Descriptions  of  each  field  (including  units  used) 
separated  by  semi-colons  (;) .  The  length  of  each  file  is  approximately  1840 
lines . 


Column : 

Format : 

Description : 

Units : 

[1-2] 

21 

Two-digit  Year 

years 

[4-6] 

31 

Day  of  the  Year 

days 

[8-15] 

8.3F 

Second  of  the  Day  GPS  time 

sec 

[17-19] 

31 

Latitude  Center  of 

Current  Averaging  Bin 

deg 

[21-29] 

9.5F 

Geodetic  Latitude 

deg 

[31-40] 

10. 5F 

Longitude 

deg 

[42-48] 

7.3F 

Satellite  Height 

km 

[50-56] 

7.4F 

Satellite  Local  Time 

hours 

[58-66] 

9.5F 

Quasi-Dipole  Latitude 

deg 

[68-77] 

10. 5F 

Magnetic  Longitude 

deg 

[79-85] 

7.4F 

Magnetic  Local  Time 

hours 

[87-98] 

12 . 6E 

Neutral  Density 

kg/mA3 

[100-111] 

12 . 6E 

Density  @  400km 

kg/mA3 

[113-124] 

12 . 6E 

Density  @  410km 

kg/mA3 

[126-137] 

12 . 6E 

NRL-MSIS  Density  at 

Satellite  Height 

kg/mA3 

[139-150] 

12 . 6E 

Density  Uncertainty 

kg/mA3 

[152-153] 

21 

Number  of  Data  Points 
in  Current  Averaging  Bin 

[155-156] 

21 

Number  of  Points  that 

Require  Interpolation 

[158-162] 

5.3F 

Coefficient  of  Drag 

Averaged  over  Bin 

%%%%  Matlab  and  Netcdf  Format: 

Matlab  and  Netcdf  formats  are  identical.  The  variables  below  are 
structure  arrays  containing  '.data',  '.units',  and  '.long  name'.  For  most  of 
the  variables,  the  size  of  the  ' .data'  structure  is  n  by  1,  where  n  is 
approximately  1840. 
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Variable  Name: 

Size : 

.long  name: 

. units : 

Year . data 

lxl 

Two-digit  Year 

years 

Doy . data 

lxl 

Day  of  the  Year 

days 

Sec . data 

nxl 

Second  of  the  Day  GPS  time. 

sec 

LatBin . data 

nxl 

Latitude  Center  of 

Current  Averaging  Bin 

deg 

Lat . data 

nxl 

Geodetic  Latitude 

deg 

Lon . data 

nxl 

Longitude 

deg 

Height . data 

nxl 

Satellite  Height 

km 

LocTim. data 

nxl 

Satellite  Local  Time 

hours 

Mlat . data 

nxl 

Quasi-Dipole  Latitude 

deg 

Mlon . data 

nxl 

Magnetic  Longitude 

deg 

Mlt . data 

nxl 

Magnetic  Local  Time 

hours 

Density . data 

nxl 

Neutral  Density 

3* 

> 

OJ 

D400 . data 

nxl 

Density  @  400km 

kg/mA3 

D410 .data 

nxl 

Density  @  410km 

kg/mA3 

Dmsis . data 

nxl 

NRL-MSIS  Density  at 

Satellite  Position 

(Q 

3* 

> 

OJ 

U  wy.data 

nxl 

Uncertainty 

3* 

> 

00 

Num. data 

nxl 

Number  of  Data  Points 
in  Current  Averaging  Bin 

Numlnterp . data 

nxl 

Number  of  Points  that 

Require  Interpolation 

Cd . data 

nxl 

Coefficient  of  Drag 

Averaged  over  Bin 

9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'S-9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9' 

ooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 

These  data  are  process  and  made  publicly-accessible  as  a  service  to  the 
scientific  community.  Individuals  planning  scientific  investigations  with 
these  data  should  contact  Dr.  Eric  Sutton  at  "AFRL/RVBXI  Org  Mailbox 
<af rl . rvbxiorgmailbox@kirtland . af ,mil>"  to  verify  the  status  of  possible  new 
versions  of  the  data,  and  to  arrange  coordination  and/or  collaboration  with 
our  ongoing  scientific  activities,  as  appropriate. 

This  work  is  supported  by  the  Orbital  Drag  Environment  (ODE)  Program  through 
the  Air  Force  Research  Laboratory  (AFRL) . 
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